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Abstract 

Peripheral collisions give access to a set of events where hot quasi- 
projectile (QP) sources are produced exploring a large range of excita- 
tion energy. In this range evaporation and multifragmentation are both 
observed for a similar charge/size of the sources. In this work, this ex- 
perimental fact will be described in terms of coexistence of two classes 
of events in the first order phase transition formalism. For this, a de- 
tailed study of the experimental correlation between the size/charge 
of the biggest cluster/fragment (Zi) and the excitation energy (E*) 
is made. Making the parallel with the same correlation derived from 
the canonical ensemble description, a first value of the latent heat and 
boundaries of spinodal and coexistence zone are extracted. 



1 Introduction 



The distribution of an order parameter of a finite system which passes 
through a first order phase transition is expected to be bimodal [1]. For 
hot nuclei, this bimodality is an echo of two classes of events which undergo 
evaporation or multifragmentation. To be observed the populations of the 
two de-excitation channels has to be roughly equivalent. This is conditioned 
first of all by a sufficient exploration of the phase diagram. Concerning 
this point, semi-peripheral collisions of heavy ion collisions produce in the 
outgoing channel hot quasi-projectile (QP) nuclei with a large range of exci- 
tation energy. The second condition to observe bimodality is the equivalent 
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population of the two phases. In heavy ion colhsions, impact parameter dis- 
tributions mainly govern the distribution of dissipated energy and largely 
favour low dissipation - big residue - evaporation channel class of events. In 
order to get rid of this experimental fact, previous studies p| proposed to 
sample QP events with the transverse energy of the light charged particle 
emitted on the quasi-target side. This dissipation sampling was related to 
canonical temperature sampling one can perform in canonical ensemble to 
define the transition region where bimodality is observed. By this way, they 
indeed observe a clear transition between a dominant evaporation-like decay 
mode, with the biggest cluster much heavier than the second one, and a dom- 
inant fragmentation mode, with the two heaviest fragments of similar size. 
A similar behavior has been reported in |3j. Different physical scenarios have 
been invoked to interpret the phenomenon: finite-system counterpart of the 
nuclear matter liquid-gas phase transition [21 11] , Jacobi transition of highly 
deformed systems [5J, self-organized criticality induced by nucleon-nucleon 
collisions [6j. In [3j, the two decay modes were associated to different exci- 
tation energies, suggesting a temperature-induced transition with non-zero 
latent heat. The qualitative agreement between refs. [21 [3] suggests that 
bimodality is a generic phenomenon. However, differences between the two 
data sets subsist, and trigger or selection bias cannot be excluded. To dis- 
entangle between the different scenarios, it is necessary to control the role of 
the entrance channel dynamics and establish if the transition is of thermal 
character. 

For Au-|-Au collisions between 60 and 100 A.MeV performed at GSI 
and collected with the INDRA multidetector, a set of events covering an 
excitation energy range of [1;10] A.MeV with size/charge variation of the 
reconstructed source around 10% can be selected. In refs [21 E], two methods 
of selection based on the kinematical properties of events are proposed in 
order to remove from the studied set of events, so called dynamical events 
containing a contribution from the mid-rapidity region. In |9], a detailed 
study of the experimental correlation between the charge/size of the biggest 
fragment/cluster (Zi) and the excitation energy {E*) is performed for these 
sets of events. 

2 Bimodality and latent heat of Gold nuclei 

The canonical correlation of two order paremeters can be defined as 
Pcan{E*,Zi) = W{E*,Zi) e-^'^'Z^^ where W is the density of states, /3 
is the canonical temperature and Zf^ is the associated partition sum. The 
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Figure 1: top panel, from left to right: renormalized distribution of the charge of 
the biggest fragment (Zi), evolution of the mean value and the standard deviation 
of Z\ with the excitation energy (E* ); the black squares correspond to the data, 
the red line to the canonical ensemble function; bottom panel, from left to right: 
renormalized 2D experimental and canonical correlation between Z\ and E* , Chi2 
estimator obtained from the fit procedure for each cell of the correlation. 



bimodal shape of this correlation is directly connected to the residual convex 
intruder in the entropy <fS/dE'^ > with S{E*,Zi) = In W{E*,Zi) [I]. 
Studying the shape of Pcan, one can localize two generic regions on the 
phase diagram: the coexistence region correponding to the region between 
its two maxima (and directly the heat capacity which is simply their rela- 
tive distance on the E* direction) , and the spinodal region where the second 
derivative of Pcan is positive {cfS/dE'^ = cfln Pcan/dE^). This theoretical 
description has to be brought closer to experimental case and to access infor- 
mation from the bimodality signal, a procedure has been proposed in ref. [1]. 
In experimental data, the distribution of excitation energy is mainly gov- 



erned by the impact parameter distribution, low dissipations are favoured. 
Measured excitation energy distributions cannot be directly related to what 
we can expect for the distribution of an order parameter. The idea is to study 
simultaneously two observables E* and Zi which are good candidates for 
order parameters of the phase transition which occurs in nuclei [iSj . Consid- 
ering that the density of states of a nuclear system exploring a large range of 
excitation energy exhibits the generic convex intruder, we can expect to un- 
derline consequent bimodality by getting rid of entrance channel effects. We 
consider the experimental correlation Pexp{E*,Zi) cx W{E* , Zi)gexp{E*) 
the simple product between the density of states and experimental weight 
which includes entrance channel effects, detection and selection bias. To 
put face to face Pexp and Pcan, a renormalization of the two correlations 
under the constraint of a flat distribution of excitation energy is made. 
It gives = Pcan,exp{E*,Zi)/Pcan,expiE*). Boltzmann factor 

and partition sum for the canonical expression and gexp for the experimental 
expression do not contribute anymore. We can focus directly on the density 
of states: W{E* , Zi)/W{E*), and compare the two normalized correlations. 
The major assumption to apply this procedure is the E* dependence only of 
the experimental weight Qexp- This means that for a given value of E* the 
associated distribution of Zi is not biased by the entrance channel and/or 
detection and selection filter. Comparing the micro canonical distribution of 
Zi for three bombarding energies, we observe no difference in the E* range 
[2,7] A.MeV where both accuracy of the calorimetry procedure to deduce 
excitation energy, detection performance and selection of events of interest 
are almost under control. This excitation energy range is chosen for the 
next step, to obtain information related to density of states of nuclei using 
fit procedure between canonical and experimental renormalized correlation. 
For this purpose, in ref. [4\, the author proposes an analytical expression 
for the canonical correlation using the double saddle point approximation. 
Considering two phases, it consists of a Taylor expansion around the cen- 
troid associated to each phase. This leads to the convolution of two gaussian 
distributions with 11 parameters : 8 mean values and standard deviations 
(2 phases, 2 dimensions), the correlation factor between Zi and E* and the 
population of the two phases. Then the renormalization is applied and the 
obtained correlation and comparison procedure can be made with data. In 
order to minimize the number of parameters, we decide to calculate the 
correlation factor from the data and to fix the respective population of 
the two phases in order to fullfill the transition point description. In [9], 
the comparison is made for two bombarding energies of Au-|-Au reaction 
at 80 and 100 A.MeV with two selection methods described in refs. [21 [7]. 



Boundaries 


Signals 


Liquid 


Gas 


Spinodal zone[lOj 


Negative Heat Capacity 


[2.0;2.5] 


[5.5;6.5] 


Spinodal zone 


Bimodality 


[1.5;3.7] 


[6.5;7.0] 


Coexistence zone[9j 


Bimodality 


[1;1.7] 


[8.5;10.4] 



Table 1: Information on the phase diagram region deduced from the study of phase 
transition signals in quasi-projectiles of Gold. 



For the four fit procedures, we obtained 4 sets of parameters and localized 
E* and Zi centroids and proposed a first estimation for the latent heat: 
AE = Eg-Ei = 8.1{±0A)stat.{+l.2 - 0.9) 

syst.- Employing, respectively, 
two different selections (two different bombarding energies) allows to put 
systematic (statistical) error bars on these results. From the general point 
of view, this exhaustive study confirms the robustness and generic properties 
of the bimodality signal in quasi-projectile fragmentation data. On fig. [H 
best reproduction of the correlation of one data set is shown. 

To go further, a complementary information can be derived from these 
analyses. From the raw Pcan{E*,Zi) = W{E*,Zi) e-^^*Z^^ and look at 
the logarithm of the projection on the E* axis {(fS/dE"^ = (fin Pcan/dE"^). 
Region where the second derivative of In Pcan{E*) is positive give us bound- 
aries of the spinodal region. For one of the four set of parameters, an example 
of this derivation is shown in fig. [2l In table [H a summary of actual status of 
experimental boundaries for coexistence and spinodal regions are shown. For 
the spinodal region we have a good coherency between boundaries deduced 
using bimodality and negative heat capacity signal [lO] . 

3 Perspectives 

In this work we make an exhaustive comparison between the 2D canonical 
distribution of two order parameters E* and Zi and the experimental one. 
This direct comparison is possible due to a renormalization procedure and 
a double saddle point approximation to put at the same level the statistical 
E* population. We extract a value for the latent heat, boundaries of the 
coexistence zone, and confirm those of the spinodal zone on the excitation 
energy axis. An important information is the evolution of these results with 
the size of the system. Varying the size/charge of the system will allow 
to investigate the influence of Coulomb forces on the stability of nuclei at 
finite temperature. It will also allow to see the minimum size of the system 
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Figure 2: fop figure, logarithme of the projection on the E* axis of the raw canonical 
2D correlation Pcan{E* , Zi); bottom figure, second derivative of the previous one. 
The dashed lines correspond to the deduce coexistence zone boundaries, the full lines 
to those of spinodal zone. 



for which the coexistence of evaporation and multifragmentation can be 
observed in terms of bimodaUty. For these purposes similar analyses on the 
Xe+Sn collisions at 65, 80 and 100 A.MeV are in progress. Another possible 
experiment could be the exclusive measurement of U+U reaction around the 
bombarding energy of 100 A.MeV. 
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